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Exotic dark matter together with the vacuum energy or cosmological constant 
seem to dominate in the Universe. An even higher density of such matter seems 
to be gravitationally trapped in our Galaxy. Thus its direct detection is central to 
particle physics and cosmology. Current fashionable supersymmetric models pro- 
vide a natural dark matter candidate which is the lightest supersymmetric particle 
(LSP). Such models combined with fairly well understood physics like the quark 
substructure of the nucleon and the nuclear structure (form factor and/or spin 
response function), permit the evaluation of the event rate for LSP-nucleus elastic 
scattering. The thus obtained event rates are, however, very low or even unde- 
tectable. So it is imperative to exploit the modulation effect, i.e. the dependence 
of the event rate on the earth's annual motion. Also it is useful to consider the 
directional rate, i.e its dependence on the direction of the recoiling nucleus. In this 
paper we study such a modulation effect both in non directional and directional 
experiments. We calculate both the differential and the total rates using both 
isothermal, symmetric as well as only axially asymmetric, and non isothermal, due 
to caustic rings, velocity distributions. We find that in the symmetric case the 
modulation amplitude is small. The same is true for the case of caustic rings. The 
inclusion of asymmetry, with a realistic enhanced velocity dispersion in the galac- 
tocentric direction, yields an enhanced modulation effect, especially in directional 
experiments. 



I. Introduction 

In recent years the consideration of exotic dark matter has become necessary in 
order to close the Universe El. Furthermore in in order to understand the large 
scale structure of the universe it has become necessary to consider matter made 
up of particles which were non-relativistic at the time of freeze out. This is the 
cold dark matter component (CDM). The COBE data suggest that CDM 
is at least 60% 0. On the other hand during the last few years evidence has 
appeared from two different teamSp the High-z Supernova Search TeamQ and 
the Supernova Cosmology Project U U which suggests that the Universe may be 
dominated by the cosmological constant A. As a matter of fact recent data the 
situation can be adequately described by a baryonic component Qb =0.1 along 
with the exotic components £Icdm — 0.3. and £!a = 0.6 (see next section for the 
definitions). In another analysis TurnerB gives Sl m = QcDftin/r^B = 0.4. Since 
the non exotic component cannot exceed 40% of the CDM El, El, there is room for 
the exotic WIMP's (Weakly Interacting Massive Particles). In fact the DAM A 
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experiment □ has claimed the observation of one signal in direct detection of 
a WIMP, which with better statistics has subsequently been interpreted as a 
modulation signal E3. 

The above developments are in line with particle physics considerations. 
Thus, in the currently favored supersymmetric (SUSY) extensions of the stan- 
dard model, the most natural WIMP candidate is the LSP, i.e. the lightest 
supersymmetric particle. In the most favored scenarios the LSP can be simply 
described as a Majorana fermion, a JipeajLjCombination of the neutral compo- 
nents of the gauginos and HiggsinosB'tilEjO. 



II. Density Versus Cosmo-logical Constant 

The evolution of the Universe is governed by the General Theory of Relativ- 
ity. The most commonly used model is that of Friedman, which utilizes the 
Robertson- Walker metric 

(ds) 2 = {dtf - R 2 {t)[ } dr) \ + r 2 {{d9) 2 + sin 2 9(d(j)) 2 )] (1) 

The resulting Einstein equations are: 

-fiV ~ 2 g ^ vR = SkGnT^v + Ag^ (2) 

where Gn is Newton's constant and A is the cosmological constant. The 
equation for the scale factor (t) becomes: 

d 2 R 4tt AttGnp A 

= -—G N {p + ?,p)R = - — (3) 

where p is the mass density. Then the energy is 

„ 1 ,dr x9 „ m, „o, A ~ k 

E = -m(— ) - G N — (4TrpR A ) H mR 1 = constant = m (4) 

2 dt R 6 2 

This can be equivalently be written as 

n k 8tt „ A 
H+ Rl = T GN > + 3 (5) 
where the quantity H is Hubble's constant defined by 

Hubble's constant is perhaps the most important parameter of cosmology. In 
fact it is not a constant but it changes with time. Its present day value is given 
by 
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H = (65 ± 15) km/s M' 1 (7) 

In other words Hq 1 = (1.50±).35) x 10 10 y, which is roughly equal to the age 
of the Universe. Astrophysicists conventionally write it as 

H = 100 h km/s M" 1 , 0.5 < h< 0.8 (8) 

Equations |||| coincide with those of the Newtonian theory with the following 
two types of forces: An attractive force decreasing in absolute value with the 
scale factor (Newton) and a repulsive force increasing with the scale factor 
(Einstein) 

mM 1 
F=-G N — 5- (Newton) , F — -AmR (Einstein) (9) 
Br 3 

Historically the cosmological constant was introduced by Einstein so that Gen- 
eral Relativity yields a stationary Universe, i.e. one which satisfies the condi- 
tions: 

Indeed for k > 0, the above equations lead to R = R c = constant provided 
that 

^AR C - y G NP R C = , ^kRl-^-G NP Rl = K (11) 

These equations have a non trivial solution provided that the density p and 
the cosmological constant A are related, i.e. 

A = 4ttG nP (12) 

The radius of the Universe then is given by 



* - < 13 > 



Define now 



Ct m — — , Qa — — , p v — „ — (" vacuum" density) (14) 

Pc Pc SnG N 

The critical density is 

Pc = 1.8xl0- 23 / l 2 ^ = 10/ l 2 ^^ (15) 
cm 6 m 6 



With these definitions Friedman's equation E — —k^ takes the form 

= (fl m + n A - l)H 2 (16) 
Thus we distinguish the following special cases: 

k > Q m + f^A > 1 O Closed curved Universe (17) 

k = TO + f^A = 1 Open Flat Universe (18) 

k < VL m + Sl\ < 1 Open Curved U niverse (19) 

In other words it is the combination of matter and "vacuum" energy, which 
determines the fate of the our Universe. 

Before concluding this section we remark that the above equations do not 
suffice to yield a solution since the density is a function of the scale factor. 
An equation of state is in addition needed, but we are not going to elaborate 
further. 

III. An Overview of Direct Detection - The Allowed SUSY Pa- 
rameter Space. 

Since this particle is expected to be very massive, m x > SOGeV, and extremely 
non relatiidstic with average kinetic energy T < lOOKeV, it can be directly de- 
tected mainly via the recoiling of a nucleus (A,Z) in the elastic scattering 
process: 

X + (A,Z) - X + (A,Z)* (20) 

(x denotes the LSP). In order to compute the event rate one needs the following 
ingredients: 

1) An effective Lagrangian at the elementary particle (quark) lejjeL ob- 
tained in the framework of supersymmetry as described , e.g., in Refs.ulij. 

2) A procedure in going from the quark to the nucleon level, i.e. a quark 
model for the nucleon. The results depend crucially on the content of the 
nucleon in quarks other than u and d. This is PjactiGularly true for the scalar 
couplings as well as the isoscalar axial coupling E3 i!£L 

3) Compute the relevant nuclear matrix elementst3 : E3 using as reliable as 
possible many body nuclear wave functions. By putting as accurate nuclear 
physics input as possible, one will be able to constrain the SUSY parameters 
as much as possible. The situation is a bit simpler in the case of the scalar 
coupling, in which case one only needs the nuclear form factor. 
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Since the obtained rates are very low, one would like to be able to exploit 
lulation of the event rates due to the earth's revolution around the sun 
l-To this end one adopts a folding procedure assuming some distribu- 
te] of velocities for the LSP. One also would like to know the directional 
rates, by observing the nucleus in a certain direction, which correlate with the 
on of the sun around the center of the galaxy and the motion of the Earth 
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The calculation of this cross section has become pretty standard. One 
starts with representative input in the restricted SUSY parameter space as 
described in the literature E3o. We will adopt a phenomenogical procedure 
taking universal soft SUSY breaking terms at Mqut, i-e., a common mass for 
all scalar fields mo, a common gaugino mass Mi/ 2 and a common trilinear 
scalar coupling Aq, which we put equal to zero (we will discuss later the influ- 
ence of nonj*ero Aq's). Our effective theory below Mqut then depends on the 



parameters 113: 

m , M 1/2 , no, a G , M GUT , h t , , h b , , h T , tan/3 , 

where a>c = 9c (#G being the GUT gauge coupling constant) and ht, hb, h T 
are respectively the top, bottom and tau Yukawa coupling constants at Mqut- 
The values of a G and Mqut are obtained as described in Refil3. For a specified 
value of tan/3 at Ms, we determine ht at Mqut by fixing the top quark mass 
at the center of its experimental range, mt(rrit) = 166GeV. The value of h T 
at Mqjjt is fixed by using the running tau lepton mass at mz, m T {mz) = 
1.746GeV. The value of hb at Mqut used is such that: 

m b{ m z)JE = 2.90 ± 0.14 GeV. 

after including the SUSY threshold correction. The SUSY parameter space is 
subject to the following constraints: 

1. ) The LSP relic abundance will satisfy the cosmological constrain: 

0.09 < fl LSP h 2 < 0.22 (21) 

2. ) The Higgs bound obtained from recent CDfEI and LEP20, i.e. m h > 113 GeV. 

3. ) We will limit ourselves to LSP-nuclcon cross sections for the scalar cou- 
pling, which gives detectable rates 

4 x 10~ 7 pb < alZiar n < 2 x 10~ 5 pb (22) 

We should remember that the event rate does not depend only on the nucleon 
cross section, but on other parameters also, mainly on the LSP mass and the 
nucleus used in target. The condition on the nucleon cross section imposes 
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severe constraints on the acceptable parameter space. In particular in our 
model it restricts tan(3 to values tan(3 — r-5p. We will not elaborate further on 
this point, since it has already appearedEJ. 



IV. Expressions for the Differential Cross Section . 



The effectivp-Lagrangian describing the LSP-nucleus cross section can be cast 
in the formLJ 

Leff = -^|{(xi7 A 75Xi) Jx + (xiXi)J} (23) 

where 

Jx = N 7x (f°+fl,T 3 + f° Al5 + f 1 Al5 T 3 )N , J = N(f° + ftr 3 )N (24) 

We have neglected the uninteresting pseudoscalar and tensor currents. 
Note that, due to the Majorana nature of the LSP, XiJ X Xi — (identically). 

With, .the above ingredients the differential cross section can be cast in the 
fonnEKMS 



du 



Mu,v) = 7TT -—^[(t s + t v —) F 2 (u) + t spm F 11 (u)} (25) 



yi spin j- > 

^spin & v o ^> spin •> Sspin f n Dill) ( I 

3(l + 4) 2 

J A 

s(u) = [Sn (o)r^ + 2^0(0)^(0)^ + nm) 2 ] (28) 

Sv = < x o Cy (29) 
(p-r(A)/ fi r {N)) 2 2 /i A-2Z, 2 .,iA), 2ri 1 277 + 1 ( 2«). 



(1 + 4)2 V A C (2 Mr &)2(l + V?) 2 (" 2 ) J 

(30) 



6 



a p x° ~ P ro * on cross-section, i = S, spin, V given by: 

a px° = a ° (fs) 2 Cm^ ) 2 (scalar) , (the isovector scalar is negligible, i.e. 

Of = CT n) 

= a 3 (/I + /A) 2 (^fl) 2 (spin) , ^ = <7 (/£ + ,^) 2 (^) 2 
(vector) 

where mjy is the nucleon mass, 77 = ra x jm,]^A, and /x r (A) is the LSP-nucleus 
reduced mass, /x r (iV) is the LSP-nucleon reduced mass and 

cr = -^(Gftojv) 2 ~ 0.77 x 1(T 38 cto 2 (31) 

Q = Q w , Qo = , 1 1? = 4.1 x 10 4 A- 4 / 3 Kel/ (32) 

where Q is the energy transfer to the nucleusr, F(u) is the nuclear form factor 
and 

F ""' {u) = g ^(oT "Mop p ' p = °' 1 (33) 

are the spin form factors (p, p are isospin indices) normalized to one at 
u = 0. Sl (^1) are t ne static isoscalar (isovector) spin matrix elements. Note 
that the quantity S(u) is essentially independent of u. So the energy transfer 
dependence is contained in the function Fn(u). Note also that S(u) depends on 
the ratio of the isoscalar to isovector axial current couplings. These individual 
couplings can vary a lot within the SUSY parameter space. 



TABLE I.: The static spin matrix elements for the light nuclei considered here. For 
comparison war-also quote the results for the medium heavy nucleus 73 Geo and the heavy 
nucleus 207 PbM 





19 p 


29 Si 


23 Na 


73 Ge 


207p b 


[Qo(0)] 2 
["i(0)] 2 

n (o)ni(o) 

f-Lexp 

Hth{spin) 1 ii exp (spin) 


2.610 
2.807 
2.707 
2.91 
2.62 
0.91 


0.207 
0.219 
-0.213 
-0.50 
-0.56 
0.99 


0.477 
0.346 
0.406 
2.22 
2.22 
0.57 


1.157 
1.005 
-1.078 


0.305 
0.231 
-0.266 
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FIG. 1.: The energy dependence of the coherent process, i.e. the square of the form factor, 
(|F(«)| 2 ), for the isotopes 19 F, 23 Na and 29 Si. The allowed range of u for the above isotopes 
is 0.011 < u < 0.17, 0.015 < u < 0.30, and 0.021 < u < 0.50 respectively. This corresponds 
to energy transfers 8.9 < Q < 140, 9.5 < Q < 190, and 9.7 < u < 230 KeV respectively. 



!_. Their ratio, however, is not changing very much. In fact actual calculations 
B show that 3.0 < 5(0) << 7.5 for 19 F, 0.03 < 5(0) < 0.2 for 29 Si and 
0.4 < 5(0) < 1.1 for 23 Na. The quantity S(u) depends very sensitively on 
nuclear physics via the static spin ME. This is exhibited in Table (jjj). As we 
can see from Table | the spin matrix elements are very accurate. This is evident 
by comparing the obtained magnetic dipole moments to experiment and noting 
that the magnetic moments, with the exception of 23 Na are dominated by the 
spin. From the same table we see that 19 F is favored from the point of view 
of the spin matrix element. This advantage may be partially lost if the LSP is 
very heavy, due to the kinematic factor fi r (A), which tends to favor a heavy 
target. The energy transfer dependence of the differential cross section for the 
coherent mode is given by the square of the form factor, i.e. |F(u)| 2 . These 
form factors for the isotopes 19 F, 23 Na and 29 Si were calculated by Divari et al 
c3 and are shown in Fig ([!]) . The energy transfer dependence of the differential 
cross section due to spin is essentially given by Fn(u). These functions for the 
isotopes 19 F 23 Na and 29 Si were calculated by Divari et alc3 and are shown in 
Fig (Q) . Note that the energy dependence of the coherent and the spin modes 
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for light systems are not very different, especially if the PCAC corrections on 
the spin response function are ignored. 




u 



FIG. 2.: The energy dependence of the spin contribution (spin response function Fn(u)) 
for the isotopes 19 F, 2S Na and 29 Si. The allowed range of energy transfers is the same as 
in Table hi In this figure we also plot Fn(u) when the PCAC effect is considered. 



V. Expressions for the Rates. 

The non-directional event rate is given by: 

dN p(0) m 

R = Rnon-dir = —7- = : a{u,v)\v\ (34) 

at m x ArriN 

Where p(0) = 0.3GeV/cm 3 is the LSP density in our vicinity and m is the 
detector mass The differential non-directional rate can be written as 

p(0) Tfl 

dR = dR n0 n~d%r = -. dcr(u,v)\v\ (35) 

m x Attin 

where da(u, v) was given above. . .. . 

The directional differential rateliJp3 in the direction e is given by : 
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dRdir = — : v.eHOv.e) — da(u,v) (36) 

m x Aitln 2tt 

where H the Heaviside step function. The factor of l/2w is introduced, since 
the differential cross section of the last equation is the same with that entering 
the non-directional rate, i.e. after an integration over the azimuthal angle 
around the nuclear momentum has been performed. In other words, crudely 
speaking, l/(27r) is the suppression factor we expect in the directional rate 
compared to the usual one. The precise suppression factor depends, of course, 
on the direction of observation. In spite of their very interesting experimental 
signatures, we will not be concerned here with directional rates. The mean 
value of the non-directional event rate of Eq. ([55]) , is obtained by convoluting 
the above expressions with the LSP velocity distribution f(v, ve) with respect 
to the Earth, i.e. is given by: 



/dR\ p(0) m f da(u,v) 3 

{-r) = 1 / f(v,v E )\v\ — <Tv 

\ du I m x ArriN J du 

The above expression can be more conveniently written as 



(37) 



(^) = ^f^) (38) 
\ du I m x Am n du 



where 



(-)=/ 



M tr.. .. -,da(u,v) s 



f(v,v E ) V ' d 6 v (39) 



^) Jy du 



After performing the needed integrations over the velocity distribution, 
to first order in the Earth's velocity, and over the energy transfer u the last 
expression takes the form 

R = Rt[l + h(a, Q min )cosa)] (40) 

where a is the phase of the Earth (a = around June 2nd) and Q m in is the 
energy transfer cutoff imposed by the detectot_In the above expressions R is 
the rate obtained in the conventional approacht£l by neglecting the folding with 
the LSP velocity and the momentum transfer dependence of the differential 
cross section, i.e. by 

P(0) 



*=^^v^P* + wys.] (4i) 

where Sj, i = S,V, spin have been defined above, see Eqs @ - @. It con- 
tains all the parameters of the SUSY models. The modulation is described 
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by the parameter h . Once the rate is known and the parameters t and h, 
which depend only on the LSP mass, the nuclear form factor and the veloc- 
ity distribution the nucleon cross section can be extracted and compared to 
experiment. 

The total directional event rates can be obtained in a similar fashion by by 



integrating Eq. (36) with respect to the velocity as well as the energy transfer 
u. We find 

Rdir = -R[(£°/47r) |(1 + hi(a, Q mm )cosa)e z .e 

- h 2 (a, Qmin) cosae y .e + h 3 (a, Q mm ) sinae x .e\ (42) 

We remind that the z-axis is in the direction of the sun's motion, the y-axis is 
perpendicular to the plane of the galaxy and the x-axis is in the galactocentric 
direction. The effect of folding with LSP velocity on the total rate is taken 
into account via the quantity t°, which depends on the LSP mass. All other 
SUSY parameters have been absorbed in R. We see that the modulation of 
the directional total event rate can be described in terms of three parameters 
hi, 1=1,2,3. In the special case of A = we essentially have one parameter, 
namely h\, since then we have hi = 0.117 and h^ = 0.135. Given the functions 
hi(a,Q m i n ) one can plot the the expression in Eq. (|4^ ) as a function of the 
phase of the earth a. 

VI. The Scalar Contribution- The Role of the Heavy Quarks 



The coherent scattering can be mediated via the the neutral intermediate Higgs 
particles (h and H), which survive as physical particles. It can also be mediated 
via s-quarks, via the mixing of the isodoublet and isosinlet s-quarks of the same 
charge. In our model we find that the Higgs contribution becomes dominant 
and, as a matter of fact the heavy Higgs H is more important (the Higgs 
particle A couples in a pseudoscalar way, whichxLoes not lead to coherence). It 
is well known that all quark flavors contribute E3, since the relevant couplings 
are proportional to the quark masses. One encounters in the nucleon not only 
the usual sea quarks (uu, dd and ss) but the heavier quarks cJx t which couple 
to the nucleon via two gluon exchange, see e.g. Drees et al 113 and references 
therein. 

As a result one obtains an effective scalar Higgs-nucleon coupling by using 
effective quark masses as follows 

m u -> f u m N , m d -> f d m N . m s -> f s m N 
rriQ —>■ /qTOjv, (heavy quarks c,b,t) 
where is the nucleon mass. The isovector contribution is now negligible. 
The parameters f q , q = u,d,s can be obtained by chiral symmetry breaking 
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terms in^-elation to phase shift and dispersion analysis. Following Cheng and 
Cheng E2l we obtain: 

f u = 0.021, f d = 0.037, f s = 0.140 (model B) 



A, = 0.023, f d = 0.034, /,= 0.400 (model C) 

We see that in both models the s-quark is dominant. Then to leading order 
via quark loops and gluon exchange with the nucleon one finds: 

/q = 2/27(1 — 

This yields: 

f Q = 0.060 (model B), f Q = 0.040 (model C) 

There is_a correction to the above parameters coming from loocs involving 
s-quarkslij and due to QCD effects. Thus for large tanp we find Ex 

f c = 0.060 x 1.068 = 0.064, f t = 0.060 x 2.048 = 0.123, f b = 0.060 x 1.174 = 

0.070 (model B) 



f c = 0.040 x 1.068 = 0.043, f t = 0.040 x 2.048 = 0.082, f b = 0.040 x 1.174 

0.047 (model B) 

For a more detailed discussion we refer the reader to Refsll 



VII. Results and Discussion 



The three basic ingredients of our calculation were the input SUSY parameters 
(see sect. 1), a quark model for the nucleon (see sect. 3) and the velocity 
distribution combined with the structure of the nuclei involved (see sect. 2). 
we will focus our attention on the coherent scattering and present results for the 
popular target 127 1. We have utilized two nucleon models indicated by B and 
C which take into account the presence of heavy quarks in the nucleon. We also 
considered energy cut offs imposed by the detector, by considering two typical 
cases Qmin = 10, 20 KeV. The thus obtained results for the unmodulated total 
non directional event rates Rt in thepcase of the symmetric isothermal model 
for a typical SUSY parameter choiceEd are shown in Fig. ra. 
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FIG. 3.: The Total detection rate per (kg — target)yr vs the LSP mass in GeV for a typical 
solution in our parameter space in the case of 127 1 corresponding to model B (thick line) 
and Model C (fine line). For the definitions see text. 



Special attention was paid to the the directional rate and its modulation 
due to the annual motion of the earth in the case of isothermal models— The 
case of non isothermal models, e.g. caustic rings, is more complicated E3 and 
it will not be further discussed here. As expected, the parameter to, which 
contains the effect of the nuclear form factor and the LSP velocity dependence, 
decreases as the reduced mass increases. 

We will focus to the discussion of the directional rates described in terms 
of to and hi, i = 1,2,3 (see Eq. ([l2])) and limit ourselves to directions of obser- 
vation close to the coordinate axes. As expected, the parameter to, decreases 
as the reduced mass increases. The quantity t° is shown in Fig. (0), for three 
values of the detector energy cutoff , Q m in = 0, 10 and 20 KeV. Similarly we 
show the quantity hi in Fig. (|^) . The quanities hi and h% are shown in Fig (^|) 
for A = 1. For A = they are not shown, since they are essentially constant 
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and equal to 0.117 and 0.135 respectively. 




A = A = 1 



FIG. 4.: The dependence of the quantity t° on the LSP mass for the symmetric case (A = 0) 
as well as for the maximum axial asymmetry (A = 1) in the case of the target 127 1. For 
orientation purposes three detection cutoff energies are exhibited, Q m i n = (thick solid 
line),Q m i n = 5 keV (thin solid line) and Q m i„ = 10 keV (dahed line). As expected t° 
decreases as the cutoff energy increases. 




A = A = 1 

FIG. 5.: The same as in the previous figure for the modulation amplitude hi. 
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h 2 (for A = 1) 



h 3 (for A = 1) 



FIG. 6.: The same as in the previous figure for the modulation amplitudes /12 and /13 for 
A = 1. 

As expected, the parameter to, decreases as the reduced mass increases. 
It also decreasres as the cutoff energy Q m in increases. We notice that t° is 
affected little by the presence of asymmetry. On the other hand hi , h 2 and h 3 
substantially increase in the presence of asymmetry. Sometimes they increase 
as the cutoff energy increases (at the expense, of course, of the total number 




iunts. For the differential rate the reader is referred to our previous work 



VIII. Conclusions 

In the present paper we have discussed the parameters, which describe the 
event rates for direct detection of SUSY dark matter. Only in a small segmant 
of the allowed parameter space the rates are above the present experimental 
goals. We thus looked for characteristic experimental signatures for back- 
ground reduction, i.e. a) Correlation of the event rates with the motion of the 
Earth (modulation effect) and b) the directional rates (their correlation both 
with the velocity of the sun and that of the Earth.) 

A typical graph for the total unmodulated rate is shown Fig. ||. We will 
concentrate here on the directional rates, described in terms of the parameters 
to,hi,h 2 and /i 3 . For simplicity these parameters are given in Figs (|J)-(||) 
for directions of observation close to the three axes x,y,z. We see that the 
unmodulated rate scales by the cos8 s , with 9 S (the angle between the direction 
of observation and the velocity of the sun). The reduction factor, f re d — 
io/(47r to) = k/(2it), of the total directional rate, along the sun's direction 
of motion, compared to the total non directional rate depends on-the nuclear 
parameters, the reduced mass and the asymmetry parameter A E3. We find 
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that k is around 0.6 (no asymmetry) and around 0.7 (maximum asymmetry, 
A = 1.0), i.e. not very different from the naively expected f re d — l/(2iv), i.e. 
k = 1. The modulation of the directional rate increases with the asymmetry 
parameter A and it also depends of the direction of observation. For Q m in = 
it can reach values up to 23%. Values up to 35% are possible for large values 
of Qmin, but they occur at the expense of the total number of counts. 

This work was supported by the European Union under the contracts 
RTN No HPRN-CT-2000-00148 and TMR No. ERBFMRX-CT96-0090 and 
HENEA 95 of the Greek Secretariat for Research. 
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